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Many radiation transport problems can best be solved by using both discrete 
ordinates and Monte Carlo techniques with a coupling between the two techniques 
occurring at a geometry interface. A general discussion of two possible coupling 
schemes is given. The calculation of the reactor radiation scattered . from a 
docked service and command module is used as an example of coupling discrete 
ordinates (DOT) and Monte Carlo (MORSE) calculations. 


Many SNAP shielding problems involve determin- 
ing neutron and gamma-ray transport in three-dimen- 
sional geometries. The solution of these problems 
by the application of Monte Carlo techniques alone 
frequently requires considerable computer time; 
however, many parts of the transport problem are 
only one- or two-dimensional and can easily be 
solved with discrete ordinates techniques. A con- 
venient approach to solving the three-dimensional 
problem is to combine the discrete ordinates and 
Monte Carlo techniques, using Monte Carlo only for 
the transport that is three dimensional. 

In coupling discrete ordinates with Monte Car- 
lo, there is generally some difficulty associated 
with the compatibility of cross-section representa- 
tion. This is due to the use of multigroup cross 
sections in the discrete ordinates codes and point 
cross sections in the Monte Carlo calculations. 
Additional problems in performing coupled neutron- 
gamma-ray calculations in the Monte Carlo part also 
arise in some cases. The use of the multigroup 
Monte Carlo code MORSE (ref. 1) alleviates both of 
these problems in that the same multigroup cross 
sections can be used in both the discrete ordinates 
and the Monte Carlo calculations, and the coupled 
neutron- gamma -ray transport can be solved as one 
problem. 

In using the results of discrete ordinates cal- 
culations as a source term for a Monte Carlo calcu- 
lation, it is desirable to use as much as possible 
of the detail that is available, or else determine 
the effect of ignoring the detail. In other words, 
the assumption of separation of energy, angular, 
and spatial variation of the radiation should be 
avoided . 
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In general, the source term for the Monte Carlo 
calculation may be distributed over a volume or may 
be a boundary or surface source. Both types of 
coupling will be discussed, and then an example will 
be given to illustrate a specific application. 

The first case to be discussed is the coupling 
of results of discrete ordinates calculations [in 
this case, DOT (ref# 2) results] to Monte Carlo 
with a volume or region being unique to both calcu- 
lations. The results of a DOT calculation of the 
fluence $(r ,z ,E,0 , *) , a five-dimensional array de^ 
fined for r,z within the volume of interest, is 
assumed to exist. In general, one would like to 
select the source term for subsequent calculations 
from a conditional probability distribution for 
particles leaving a collision as a function of each 
of the five variables. Because of the large array 
of numbers , approximations have been made in the 
past in which some of the variables are selected 
from non-conditional probabilities that are volume 
integrated. For example, the energy may be sel- 
ected from the spectrum integrated over the volume. 
Likewise, the angular distribution may be assumed 
to be the same for all points in the volume, or for 
all energies. However, complete detail can be used 
with increased data handling and storage require- 
ments . 

The first step in coupling over a spatial reg- 
ion is to convert the fluence $ to x» the density 
of particles leaving a collision. This is accom- 
plished by multiplying the fluence by the group-t©- 
group cross section: 

x(r , z,E ,6 ,<t>) = l l $(r,z,E’,e',4>) 
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T = I 2irrAr [ Az J J J wt(A,m) x (r,z,E ,0,*) 
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where x(r , z ,E g ,9 J) 2irrArAz is defined as the densi- 
ty of particles leaving a collision at r,z with 
energy E in group g per unit polar angle 0 and azi- 
muthal angle <J>. 

wt(£) = £ wt(£,m) is the solid angle segment 
m 

for each polar angle 0 , and 
wt (&,m) = the solid angle of each azimuthal angle 
4> for each polar angle index Z, 

In order to select from this distribution, one 
must form the five following distributions: 
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case, similar conditional probabilities must also 
be formed, with some minor exceptions. Consider 
the coupling surface to be the outer boundary of a 
cylinder . 

The leakage fluence, $(r ,Z ,E g , 0 ,<{>) , at the top, 
bottom, and curved surface of a cylinder is deter- 
mined as a function of energy, polar, and azimuthal 
angles, as well as the spatial variation on the sur- 
faces. The leakage fluence on these surfaces desig- 
nated as R q , Z^, and must be converted to current 
for the selection of the source term of the Monte 
Carlo calculation. Thus, there are three spatial 
terms - the radial variation at the top, the radial 
variation at the bottom, and the height variation 
at the outer surface of the cylinder. Define T, B, 
and S for the integrated leakage for the top, bottom, 
and side of the cylinder. 

T = I 2irrAr £ £ cos(0 )£ wt(A,m) 
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From these conditional probability distributions 
one may select values for the source coordinates. 

With the assumption of linear variation between 
spatial points and angular directions the discrete 
values are not carried forth in the Monte Carlo 
calculation, but instead a continuous distribution 
of coordinates is selected. Care must be exercised 
in utilizing this type of coupling so that the 
transport in the overlap region is not included 
twice. This may be accomplished using a pure ab- 
sorber in the overlap region in the Monte Carlo cal- 
culation • 


wt(£ ,m)cos<J> <KR o ,Z,E g ,0,<J>) 


with total leakage TL = T + B + S. 

each of the above terms , a four— dimensional 
conditional probability must be determined. For 
example, consider the bottom surface leakage: 
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Perhaps a more useful coupling technique is 
the coupling at a boundary of the system. For this 


675 



*SUM . 




SjS W!« M 


b x g 


wate 


wate* 


P^) 

|P(Xj)| 


Ag(<|> J 0 »E^,r ) 


x $(r ,Zg,E^ ,0 , (Jt) » 

2TrrAr|cose|wtg > m) 

BR(r)E B (E g |r)P B (0|E g ,r) 


x $(r >Zg*E g ,0 ,$) . 


The results for the top surface leakage and side 
leakage are similar . 

With probability B/TL and T/TL the source 
would be selected from the bottom and top surface, 
respectively, and with probability S/TL the source 
would be selected from the side of the cylinder. 
Then, given the spatial coordinates, an energy 
group is chosen from Eg, and then the direction 
cosines are chosen for P fi and A fi . 

The selection from the probability distribution 
given above is sometimes hindered by negative fluxes 
in the DOT-calculated source term. If there are no 
negative values, then the selection of a value x is 
selected from the distribution by choosing a ran- 
dom number R and finding j such that 

J-l i 

l P(x.) < R £ l P(V ; 

k=0 k k=0 


then 
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However, if some values of the distribution 
P(x^) are negative, then an alternate scheme is 
required. The random number is scaled by the sum 
of the absolute values of the probabilities and a 
selection is made. A weight correction is required 
which can result in negative particle weights. 

These negative weights. in general will not severely 
affect the solution, but they do, of course, affect 
the variance of the result. 


j-1 j 

l |P(X )|<R*SUM < l |P(x,)| , 
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where 

Nk 

SUM =* l |P(x, )| 
k=0 


If there is a large probability of selecting a neg- 
ative weight, then the DOT source term should be 
recalculated. 

To illustrate DOT -MORSE coupling, we have con- 
sidered a problem in which the coupling is at a 
boundary surface. More specifically, the problem 
consists of calculating the radiation scattered 
from a docked command and service module due to the 
leakage radiation from a 600-kWt ZrH reactor (with 
its associated shielding) . The problem is illus- 
trated in figure 1. The problem which consists of 
determining the scattered neutron and gamma -ray 
dose at a detector plane due to neutrons and gamma 
rays leaking from the reactor power assembly has 
previously been investigated at Atomics Internation- 
al (refs. 3 and 4). A drawing of the reactor and 
shield assembly is shown in figure 2. 


OflNl-OWG 7< - 1,39 



FIGURE 1. -Schematic of Reactor Power System and 
Docked Service and Command Module. 


and 
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FIGURE 2. -Reactor Shield Assembly. 

There are several approaches that could be used 
in determining the dose due to radiation scattering 
from the docked module. For example, one may use 
the boundary leakage flux obtained in a discrete 
ordinates calculation of the reactor and shield 
assembly as a source term for the Monte Carlo calcu- 
lation, transport the radiation to the module, 
follow the transport in the service module, and 
estimate the dose at the detector plane. A second 
approach is to perform an adjoint calculation of the 
importance of radiation leaking from the core-shield 
assembly by Monte Carlo and then integrate the impor- 
tance function over the source of particles leaking 
from the reactor assembly which was obtained in a 
discrete ordinates calculation. Another approach 
is to determine the spatial and energy distribution 
of the fluence at a surface surrounding the module 
and use this surface as the boundary for coupling 
the discrete ordinates results to a Monte Carlo cal- 
culation of the scattering in the service module. 

This last approach will be discussed. 

A two-dimensional calculation of the reactor 
and shield assembly, as shown in figure 2, was made 
with the DOT (ref. 2) code. An Sg quadrature was 
used to determine the energy spectra in 21 neutron 
and 18 gamma-ray groups. The radiation field at 15 
radial locations along a plane at the top of the 
service module was obtained with SPACETRAN (ref. 2) 
using the leakage boundary fluence. SPACETRAN per- 
forms a ray-tracing calculation from each spatial 


mesh point to the detector point. The intensity 
and energy distribution of both neutrons and gamma 
rays was determined for each of these 15 locations. 
See Table I for energy group structure and fluence- 
to-dose conversion factors. The angular distribu- 
tion of the source was determined by assuming that 
all the radiation was coming from the center of the 
reactor core located approximately 90 feet away. 

(The problem is not expected to be very sensitive 
to the small angular variation of the radiation 
incident on the module.) This radiation field as 
a function of space and energy was used as the source 
term for Monte Carlo. 


Table I. Energy Group Structure and 
Fluence-To-Dose Conversion Factors 


Group 

Upper 

Energy (eV) 

Fluence-to-Dose 
Conversion Factors 
(mr em/hr)/ 
(part/cm* /sec) 

Neutrons: 

1 

1.4918 (+7)* 

1 . 5000 (-1) 

2 

1.0000 (+7) 

1.5000(-1) 

3 

6. 7G32(+6) 

1.3700(-1) 

4 

4 . 4933(+6) 

1 - 3 200 (—1 ) 

5 

3. 0119 (+6) 

1 . 3100 (-1) 

6 

2. 0190 (+6) 

1 . 2500(-l) 

7 

1.3534 (+6) 

1 . 1600(-1) 

8 

9.0718(+5) 

1 .0600 (-1) 

9 

5.5023(+5) 

7. 5 700 (-2) 

10 

3. 3373(+5) 

5. 51 00 (-2) 

11 

2. 0242 (+5) 

4.0100(-2) 

12 

1.2277 (+5) 

2.4500(-2) 

13 

4. 0867 (+4) 

8. 5000 (-3) 

14 

1.1709 (+4) 

5. 0000 (-3) 

15 

3. 3546(+3) 

5.0000(-3) 

16 

7.4852(+2) 

5. 0000 (-3) 

17 

1.67 02 (+2) 

5.0000<-3) 

18 

3. 7266 (+1) 

5. 0000 (-3) 

19 

8.3153(0) 

5.0000 (-3) 

20 

1.8554(0) 

5. 0000(-3) 

21 

4. 1399 (-1) 

3.7500(-3) 

Gamma Rays : 

22 

1.0000 (+7) 

9.8000(-3) 

23 

8. 0000 (+6) 

8. 5000 (-3) 

24 

7. 0000 (+6) 

7. 6000 (-3) 

25 

6. 0000 (+6) 

6.7000(-3) 

26 

5. 0000 (+6) 

5. 8000 (-3) 

27 

4. 0000 (+6) 

5. 0000 (-3) 

28 

3. 5000 (+6) 

4. 5000 (-3) 

29 

3. 0000 (+6) 

4. 0000 ( -Ti 

30 

2. 5000 (+6) 

3. 5000 (-3) 

31 

2.0000 (+6) 

3.0000(-3) 

32 

1.6000 (+6) 

2. 4000 (-3) 

33 

1.2000 (+6) 

2. 0000 (-3) 

34 

9. 0000 (+5) 

1 .5000 (-3) 

35 

6. 0000 (+5) 

1.0500 (-3) 

36 

4 . 0000 (+ 5 ) 

6 . 0000(-4) 

37 

2. 10 00 (+5) 

2. 8000 (-4) 

38 

1 . 2000(+5) 

1.4000(-4) 

39 

* 

7. 0000 (+4) 

4. 0000 (-4) 


Read as 1.4918 x 10 7 . 


At a plane just above the module, the intensity 
variation was determined to be that shown in figure 
3. The energy distribution, or the number in each 
of the 21 neutron groups and 18 gamma-ray groups , 
is shown as a cumulative distribution for three 
different radii in figure 4. These curves illus- 
trate the distributions used for the selection of 
the source in the MORSE Monte Carlo calculation. 

By using MORSE, it was possible to determine 
the contribution to the scattered dose of the 
incident neutrons and gamma rays as well as the 
secondary gamma rays produced in the module. 
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The geometry employed in the calculation is 
illustrated in figure 5, which shows boundary cross- 
ing events in the transport of neutrons in the mod- 
ule. The radial and longitudinal location of the 
crossing events is shown. This collision density 
plot is an output from MORSE and gives the boundary 
crossings for approximately 1000 source particles. 

The typed numbers give the density of the homogen- 
ized parts of the service module, and Table II 
gives the material composition for each region. 
Figure 6 shows the density of collisions of neutrons 
in one particular problem. The rotational symmetry 
of the module is taken into account in the plot, 
and there are few collisions in the center. The use 
of this type of on-line plot has been useful in 
determining where the events are occurring and the 
effects of importance sampling. 


FIGURE 3. -Spatial Distribution of Radiation at 

Source Plane Adjacent to the Service and Command 


Module . 


■ *U K../CC 
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FIGURE 5. -Boundary Crossing Events From a MORSE 
Calculation. The density of the three material 
regions is given. 


Table II. Material Composition for the 
Service and Command Module 



Weight 

(lbs) 

Density 

(gm/cc) 


Composition 

Cwt %) 


Medium 

H 

0 

A1 

Ti 

Fe 

1 

13,500 

0.612 

3.3 

46.7 

20.0 

17.0 

. 13.0 

2 

20,800 

0.193 

3.3 

46.7 

20.0 

17.0 

13.0 

3 

200 

1.36 





100.0 


FIGURE A. -Cumulative Probability for Energy 
Distribution of Neutrons and Gamma Rays at 
Three Locations Along Source Plane. 
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FIGURE 6. -Collision Density Plot for Neutrons 
Scattering in the Service and Command Module. 
Numbers 1-9 include the number of collisions in 
the intervals with 0 indicating 10 or more 
collisions. 


The contribution of scattered neutrons and 
gamma rays to the dose at the detector plane (117 
ft from the reactor) is shown in Table III (about 
5% of the gamma-ray dose is due to secondary gamma 
rays produced in the service module). The total 
scattered dose of 0.34 mrem/hr should be compared 
to the direct dose of 0.70 mrem/hr (0.66 mrem/hr 
gamma rays and 0.04 mrem/hr neutrons). These results 
are compared with values calculated at Atomics 
International for various locations on the detector 
plane. There is surprisingly good agreement be- 
tween the Monte Carlo and single-scattering results. 


Table III. Scattered Dose at the Detector Plane 
(mrem/hr) 


Detector 
Location (ft) 

Neutron Dose 

Gamma-Ray Dose 

1 Total 







Along 

Axis 

1 Axis 

Single 

Scattering* 

Monte 

Carlo 

Single 

Scattering 

Monte 

Carlo 

Single 

Scattering 

Monte 

Carlo 

0 

0 

0.228 

0.29 

0.063 

0.05 

0.30 

0.34 

11 

0 

0.73 

0.31 

0.27 

0.08 

1.0 

0.39 

0 

11 

0.27 

0.39 

0.09 

0.04 

0.36 

0.43 


* 

With self-absorption. 
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As a check on the adequacy of the SPACETRAN cal- 
culated spatial distribution at the service module 
location, a DOT calculation with a 100-angle quad- 
rature biased in the downward direction was made. 
Figures 7 and 8 show the spatial distribution at 
the source plane and the energy distributions. 

There are fairly large differences in these results 
compared to those in figures 3 and 4 ; however, the 
scattered dose contribution does not change more 
than 15%. 



RADIUS tCMJ 


FIGURE 7. -Spatial Distribution of Radiation at 
Source Plane Adjacent to the Service and 
Command Module Calculated With a 100-angle 
Asymmetric Quadrature. 



FIGURE 8. -Cumulative Probability for Energy 
Distribution of Neutrons and Gamma Rays at 
Three Locations Along Source Plane Calculated 
With a 100-angle Asymmetric Quadrature. 
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